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Abstract 

Astrophysical bounds on the cosmological constant are examined for 
spherically symmetric bodies. Similar limits emerge from hydrostatical 
and gravitational equilibrium and the validity of the Newtonian limit. 
The methods in use seem to be disjoint from the basic principles, however 
they have the same implication regarding the upper bounds. Therefore we 
will compare different inequalities and comment on the possible relation- 
ship between them. These inequalities are of relevance for the so-called 
coincidence problem and for the bound of the cosmological constant which 
comes surprisingly close to the 'experimental' value. 
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1 Introduction 



The cosmological constant A has a history as long as General Relativity itself: 
introduced to make the universe static, abandoned after the discovery of the 
Hubble expansion and introduced recently to explain the acceleration of the 
Universe [1,2]. It is true that the latter phenomenon can, in principle, be also 
accounted for by other explanations like quintessence [3,4], Chaplygin Gas [5], 
mixture models [6-9] and modified gravity [10] etc, but the simplest explanation 
for an acceleration can be obtained from a positive cosmological constant with 

A = 8irGp vlic , (1) 

Pvac ~ (0.7 - 0.8)p crit . 

This value is surprisingly close to a prediction made long ago by Zel'dovich [11- 
14], namely, 

A w G 2 m 6 , (2) 

where m is a typical hadronic mass scale (0.15 — 1 GeV). Indeed, combining 
this with H = m 3 [15,16], one gets the result that the Hubble parameter is 
dominated by the cosmological constant A. 

Cosmological observations [17-19] give strong indications of the presence of a 
positive cosmological constant which would mean that the universe is of de Sitter 
type. Moreover, recent studies [20] of X-Rays hint toward a constant density 
of the dark energy which is thought responsible for the acceleration, also in 
agreement with a positive A. The use of the different distance measurements in 
a ACDM model (models of cold dark matter with A) also shows strong evidence 
for a positive cosmological constant (0.47 < flA < 0.82, Qa = Pvac/Pcrit) [21]. 

Such a relatively large cosmological constant has astrophysical implications. 
This was realized as early as 1939 by Gamow and Teller [22]. The situation 
is best paraphrased in [23] which summarizes the work of several authors in 
the 30's and 40's: "The essential difficulty with a relativistic theory in which A 
[the Cosmological Constant] is positive is that of accounting for the formation 
and condensation in terms of gravitational instability; for, to use the 'force' 
metaphor, the present expansion indicates that the force of cosmic repulsion 
exceeds those of gravitational attraction [22]. This is not likely to disturb the 
stability of systems (such as the galaxy) of high average density, but it is likely 
to prevent new condensation in regions of low density." 

Indeed, several studies have corroborated this conclusion and many well 
known systems are reconsidered with the inclusion of the formerly omitted A- 
term. Another, more mathematical reason to include A is the richer asymptotic 
structure of these spacetimes and moreover the AdS/CFT correspondence [24] 
motivating much research. 

Examples of recent research are the Einstein- Yang-Mills- A system [25] , (n + 
4) dimensional EYM-A theory [26], perfect fluid solutions with cosmological 
constant [27-30], consistency with the Newtonian limit [31,32], gravitational 



2 



equilibrium via the virial theorem [33] and also Einstcin-Cartan and Einstein- 
Cartan-Dirac theory with cosmological constant [34, 35] . It is remarkable that 
these examples have something in common, namely, in each of them there exists 
an upper bound of the cosmological constant. 

Gravitational equilibrium and perfect fluid considerations surprisingly pre- 
dict the same upper bound in terms of the mean density of an astrophysical 
object, A < 4irGp. Hence it seems that there is some relation between the 
existence of exact solutions of Einstein's field equations and the virial theorem. 
In the EYM- A system, one finds A < A crit . In the (n + 4) dimensional EYM- 
A theory it was found that the cosmological constant is bounded, such that 
A(n+4) < 5 2 /2G(„ + 4) and constant density solutions exist for A < As, for larger 
values they are unphysical. Finally we mention that supergravity [36] even for- 
bids a positive value and therefore predicts A < 0. It should be noted that a 
lower bound never occurs. 

In [33] it was found that the more the shape of the object deviates from 
spherical symmetry, the more difficult it is to reach the gravitational equilibrium 
in presence of a positive cosmological constant for low density astrophysical 
objects. This result was based on the assumption of gravitational equilibrium. 
Of course, equally valid and useful concept is the hydrodynamical equilibrium 
combined with Einstein's equations. These two concepts of equilibrium give us 
information on A in terms of inequalities. A third source of useful information 
on A is the Newtonian limit, which also gives us an inequality. 

In this paper we will compare different inequalities and point out some inter- 
esting coincidences. We will use these results to infer on certain static properties 
of astrophysical objects. Unless otherwise stated we put G = c = h = 1. 

2 Newtonian hydrostatic equilibrium 

In order to be self consistent, some basic relations are reviewed as shortly as 
possible. 

In Newtonian astrophysics a spherically symmetric object is stable if the 
gravitational and matter pressure are in equilibrium. This condition leads to 



which is sometimes called [37, p. 301] the "fundamental equation of Newtonian 
astrophysics" . The mass function is as usual defined by 




(3) 




(4) 



Furthermore let the mean density up to r be defined by 



p(r) = 



3 m(r) 



(5) 



4ir r 3 
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which then has the following properties 



Pc = Pc, p(r) > p(r) V r , (6) 

where p c is the central density. If a density profile p(r) for the astrophysical 
object is prescribed, one can firstly use Q to find the mass function and secondly 
integrate ^ to find the pressure P(r). 

However the most physical starting point is to prescribe an equation of state 
p = p(P). First we rewrite ([3]) with the help of ([5]) and find 

P'(r) = -r^(47rp(r)-A). (7) 

Then, differentiating the mean density ([5|) leads to 

p'(r) = *(p(P(r))-p(r)). (8) 

Thus the two equations ([7]) and ([8|) form a system of differential equations in 
P(r) and p(r). 

For any physically reasonable astrophysical object, the pressure and den- 
sity must be monotonically decreasing functions of the object's radius. Hence 
negativity of P'(r) from implies 

A < 4vr ( o b , (9) 

since by virtue of (|6|) the mean density is a monotonically decreasing function 
having its minimum at the boundary of the stellar object, where it is just the 
mean density of the stellar object. The boundary is defined by the radius R for 
which the pressure vanishes, P(R) = 0. 

If, for the moment, we assume p(r) = p = const., then gives 

P'(r)=-r-|(4 7 rp-A) ! (10) 

and is, as before, well defined if and only if A < 4irp. However, if A = iirp 
we find that the pressure must be constant for all r. Hence one may say that 
this is the Newtonian analogue of the Einstein static universe. Note that the 
cosmological constant in the Newtonian case is independent of the pressure as 
is expected from general grounds since in Newton's theory of gravity pressure 
does not contribute to the energy momentum tensor. 

Moreover if the pressure increases near the center, it automatically diverges. 
This changes considerably in the general relativistic case. 



3 General relativistic hydrostatic equilibrium 

Most parts of the former section can still be applied to the general relativistic 
case. The only difference is that the differential equation for the pressure ([7]) is 
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replaced by the Tolman-Oppenheimcr-Volkoff equation [38,39] with cosmological 
constant (TOV-A) and is given as 



i,_ 3 T p(r) Af p(r)r2 - \±r* ) ' ( ' 

where the velocity of light has been included by A = 1/c 2 to facilitate the 
comparison with the Newtonian case. The Newtonian limit [40], A — > 0, of 
reproduces 0. We will use geometrical units henceforth. 

The TOV-A equation together with (jHJ form a system of differential equa- 
tions if an equation of state p = p(P) is specified. 

The right hand side of (| 1 1 1) is well defined if the denominator is 

V 2 ■= 1 - y P{r)r 2 - |r 2 > Vr, (12) 

where we introduced a new variable y. The interesting question is whether one 
can find a general relativistic analogue of the upper bound of the cosmological 
constant ©. The answer to this question is affirmative, and it is related to fact 
that the TOV-A equation be well defined. 

One can derive [27,42,43] the following inequality 

127rP(r)+47rp(r)- A 
V ~ 12nP c + 4np c -A ' (13) 

which can be used to arrive at the desired equation. Since the boundary of any 
astrophysical object is defined by P = 0, the variable y is well defined up to the 
boundary if the remaining terms of the numerator satisfy 

A < 47rp b . (14) 

Therefore, the Newtonian upper bound of the cosmological constant could ex- 
actly be reproduced in the general relativistic case. 

Before continuing, let us shortly review qualitatively how one derives (fT3|) . 
The Einstein field equations for a static and spherically symmetric perfect fluid 
are three independent differential equations for four free functions. Hence one 
can specify one of them. The most physical one is to prescribe an equation of 
state p — p(P). The TOV-A equation can be extracted from this system by 
eliminating the free function g rr if the line element is given in its usual form. 
On the other hand, one can also eliminate the pressure from the second 
field equation and the conservation equation. After introducing Buchdahl vari- 
ables [42] and assuming that the pressure is a decreasing function near the center 
and prescribing a monotonic equation of state one arrives at (|13[) . Without the 
cosmological constant one can directly read it off from [42,43]. The cosmological 
constant is incorporated by using the standard substitution p — ► p + A/8tt and 



Ms 2 = -e"Mdt a + e a Mdr 2 + r 2 {dd 2 + : 
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P — > P — A/87T which leads again to equation (fT5| . For an explicit derivation 
with the cosmological constant see [27, p.42-p.47]. 

One must be slightly careful with the above, since in case of a constant 
density distribution p — po one can integrate the TOV-A equation exactly and 
no upper bound A = Airp arises, which is expected since the Einstein static 
universe with non- vanishing pressure is also a solution of the field equations. The 
interior metric has the geometry of a three sphere which is conformally flat [44] 
and one just has to deal with a harmless coordinate singularity. See [28] for 
a complete discussion of constant density solutions with cosmological constant 
where it is also shown that in contrast to the Newtonian case one can construct 
solutions with increasing pressure near the center. 



4 Buchdahl inequalities 

There are essentially two ways of solving the field equations, as already men- 
tioned in the previous section. With the help of Buchdahl variables [42] one 
has a very convenient scheme to compare a general monotonically decreasing 
solution with a constant density solution. The latter is defined by the boundary 
of the mean density of the general solution. 

To perform the described procedure, Buchdahl [42] originally assumed the 
existence of a global solution, which was later proven to exist in [43]. Up to 
A < 4irp~b one can use the methods of [43] to prove the existence of a global 
solution with non- vanishing cosmological constant [27, Theorem 5]. 

The derivation of the Buchdahl inequality is quite involved, therefore we only 
state the result and refer to [42] and [45, p. 367] without and [27, Theorem 6] 
with the cosmological constant. The generalized Buchdahl inequality is given 

by 



- ^p b m - -R? > — . (15) 

3 H 3 " 3 12Trp b V ' 

For constant density solutions, generalized Buchdahl inequalities can be found 
in [28,46]. In [46] the surface redshift with A was derived from (fT5")) . 
For A = this inequality reduces to 



l-^p b R 2 >\, (16) 
which with M — (47r/3)pi? 3 leads to its most popular form 



I-f >i, 2„<? fl , (1 7) 

from which one reads off the well known result that radii of static perfect fluid 
spheres are larger than the Schwarzschild radius of the corresponding mass. 
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Let us for the moment assume that a global solution exists independent of 
the cosmological constant, so we assume existence for A > Airpb- From the 
generalized Buchdahl inequality (|15p one finds instead 



R 2 < 



3 \* iirp b ) 



(18) 



47T ( O b 



which after some algebra can be written as 



A 



9M 



(19) 



16-n-pb 




Since the right-hand side is bounded by 1, one simply finds 



A < l6TTp b , 



(20) 



which is in full agreement with the results from constant density solutions [28], 
where solutions exist up to a cosmological constant Ag, given by 



Therefore it seems that properties of the constant density solutions in [28] are 
also a feature of general solutions with any given equation of state. However, 
existence of solutions with A > Airpb so far could not be proved. Nonetheless we 
conjecture their existence based on the generalized Buchdahl inequality, similar 
to Buchdahl's original assumption. Moreover one may use pb = 3M/(47ri? 3 ) 
in (TT8"]) which yields 



and reduces to (fTT)) for A = 0. The above equation may also be regarded as the 
generalized Buchdahl inequality since it gives an upper bound on the mass. It is 
tempting to make a Taylor expansion with respect to A in equation (|22p , trying 
to extract some further information. However, by looking at (fT8|) and (fl9|) one 
sees that nothing can be found since both equations are linear in A. Moreover 
it should be noted that the square root appearing on the right hand side of (|22[) 
is well defined if 



is satisfied by the object's radius R. Inserting the highest possible radius (123| 
in (f22| we obtain 



A s = Wirp + 24P C . 



(21) 




(22) 




(23) 




(24) 
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It is now worthwhile to compare the inequalities (f2"5)) and with constraints 
arising from the Newtonian limit. The Newtonian approximation has, of course, 
a certain range of validity. In the static, non-relativistic limit, the requirement 
of weak fields results into two strong inequalities of the form 

» R » M , M max = |V2~4= » Af , (25) 

valid up to small corrections of the order M/M max and (M/M max ) 2 . What is 
striking, is the similarity of inequalities (f2"5| and (f2"l)) to the restrictions (|2"5"|) of 
the Newtonian limit. Disregarding different numerical factors which are of the 
order of unity and the fact that (f2"S"]) are strong inequalities, we see that both are 
essentially the same. Inequalities (|2"5)l and (fM)) can be understood as constraints 
on A to keep the object in hydrostatic equilibrium from which they are derived. 
On the other hand, ([23)) ensures that the gravitational fields are not too strong. 
Hence in both cases one expects a restriction on a positive cosmological constant. 
What is surprising, however, is the fact that these restrictions are so similar in 
both cases. We will discuss further relevance of equations (|2"5|) - (l2~5l) in section 
8. 




5 Lessons from constant density solutions 

In this section, static and spherically symmetric constant density perfect fluids 
are shortly revisited. A new inequality is derived which in the constant density 
case is shown to reduce to an equality. In the subsequent section this is used to 
extract information on the cosmological constant which is encoded in a cubic 
equation in A. 



5.1 With equation of state 

Assume than a constant density and an equation of state are prescribed, for 
instance, the equation of state of a polytrope 

P = np 1 . (26) 

In this case we have P'(r) = 0, that is, we also have a constant pressure (we can 
take this as an approximation in the sense that P'(r) « for a slowly varying 
density or pressure profile). Using this in the TOV-A equation leads to 

1 2 m / 4irp+12irP - A % 



p. (r) = - r 3 ( , + P)| - f¥ _ ? - | = 0, (27) 



and implies an equation for the constant density p: 

Amp + YhxKpi - A = . (28) 
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Note that with vanishing cosmological constant A = 0, there would not be any 
acceptable solution for a positive energy density. Thus, we can say that in the 
presence of a cosmological constant, a system with constant density and constant 
pressure can be in hydrostatical equilibrium; the Einstein static universe is 
exactly such a system. In the case with A = 0, such a system could not be 
in hydrostatical equilibrium; sooner or later the system will collapse under it's 
own gravity. From (|28|) we then can find p = p(k, 7, A). To illustrate the 
consequence of equation (|28[) let us assume that 7 = 1. It then follows that an 
astrophysical object with p ~ p crit is stable. Such stable objects would have 
the lowest possible density among stable astrophysical objects and could be e.g. 
superclusters. 

5.2 Without equation of state 

Let us assume that no equation of state is specified and consider the inequality 



The functions y(r) and ((P(r)) are the Buchdahl variables [42], already men- 
tioned above: 




(29) 




(30) 



For a constant density distribution one finds 



C(r) 



Pc + P 



(31) 



P(r) + p ' 



and equation (|2"9"|) becomes an equality 




P(r) + p V 12ttP c + 4irp c - A 



P c + p / 12nP(r) + Anp(r) - A 



) 



(32) 



Solving for the pressure and using p c = p — p we find 



P(r) = p 



" (I-/?) (P c + p[l-y(r)])-3P c y(ry 
_ 3P c (y(r)-l) + p(y(r)[l-/?]-3) . 



(33) 



where (3 = A/Anp and moreover 



Pc = p(l - 0) 



y(R) - 1 



(34) 



3y(R)-(l-0) ' 
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Finiteness of the central pressure implies that the denominator must be well 
defined, hence greater than zero. In the simplified constant density case this 
implies 

y(R)>±(l-(3), (35) 

which is just the Buchdahl inequality (fTS"|) . 

To justify why in the constant density case the inequality (j2"9"j) becomes an 
equality, one can compare with the solution of the TOV-A equation for constant 
density [28]. One finds 

where C is a constant of integration which can be fixed with the condition 
P(R) = 0. This yields 

(37) 

On the other hand, at the center r = we have 

_ (f3-l)(y(R)-l) 
Pc ~ P 3y(R)-(l~P) ' (38) 

and therefore the constant of integration becomes 

Pc+P V ' 

Solving for y(r) from (|36[) yields 

If we replace the constant C from (139(1 we obtain equation (|3"2")) and hence (f3"3")l . 
This justifies the conversion from an inequality to an equality in (|29p for constant 
density. Thus it is trivially fulfilled for constant density solutions. 



6 General solutions with equation of state 

In view of the above results, we can assume a variable density in order to find 
a more restrictive condition on A. First of all, inequality (I29|) is valid for any 
density profile; on the other hand, since it must be fulfilled at any r in order 
that the pressure be a decreasing function of the radius, we can evaluate it at 
the boundary r = R, defined by the condition P(R) = 0, which leads to 
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In this case we can safely assume an equation of state in the form of (J2SJ), so that 
the central pressure is only a function of the central density. The function £(P) 
can be integrated and becomes a function of the central density when evaluated 
at the boundary 

C(P(R)) = {fpf-^ + l) ^ = C(Pc) ■ (42) 
Equation (I41|) then becomes a cubic equation for A, namely, 

/(A) = aA 3 + bA 2 + cA + d > , (43) 
where the coefficients are given explicitly as 

a = - ^R 3 , 6=1 + 8ttP c R 2 ~ y R 2 (p ~ p c ) - ( 2 , 

c = - 2AtxP c - 48n 2 P 2 R 2 + Miz 2 pP c R 2 + 8n( 2 p - 8np c - 32n 2 P c p c R 2 

64 16 , n , , 

+ —Tr 2 R 2 p Pc - —tt 2 R 2 p 2 c , 

d = U4tt 2 P 2 - 384n 3 P 2 pR 2 - 16n 2 ( 2 p 2 + 96n 2 P cPc - 256tt 3 P c p Pc R 2 
+ 16tt 2 Pc - ^ 3 pp 2 c R 2 , (44) 

and where p = pb = p(R) is the object's mean density. Thus we have the 
following parameters R, p c , P c , pb, which can be related as follows: the central 
pressure is connected with the central density through the equation of state so 
that P c = npj. The radius of the configuration is related with the mean density 
at the surface as 

while the mass ^ is just the volume integral of the density 

/ 4irs 2 p(s)ds . (46) 
Jo 



M 



We need the density profile and the total mean density (together with an equa- 
tion of state) to solve the cubic equation (|4*3"j) for A. 

It is interesting to further exploit the inequality (|4Tj) . Although it is quite 
involved to extract information on the cosmological constant one can solve it 
for the central energy density. This yields 

Pc > Pcmin = ^T^- - 3P C + A ; (47) 

y(R) 4tt 

giving a lower bound on the central energy density. This is not surprising 
because we have upper bounds on the boundary mean density. Since the energy 
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density for astro-physical models is a decreasing function of the radius, one must 
find some lower bound, that in particular extends the boundary mean density. 

A clever choice of new variables could be of help in order to get information 
on the cosmological constant from (|4"Tj) . For example one could use the effective 
quantities mentioned in the derivation of equation (|13[) . However, the numerator 
of [|4"7]) suggests a definition of type A = p — A/ An. In future we plan to explore 
the astrophysical significance of inequality (|4ip and other versions thereof. 



7 Inequalities from virial theorem 

In what follows, the non-relativistic virial theorem is recalled and an upper 
bound on the cosmological constant is derived with surprising similarities to 
those derived before. 

Note that there exist several relativistic formulations of the virial theorem 
in literature [47-51], however, they are not complete generalizations since some 
restrictive assumptions are always needed in their formulation. 

The standard non-relativistic virial theorem reads 



d 2 I 



jk 



dt 2 



= AK jk + 2W jk + ■ • ■ , (48) 



where Ij k is the inertial tensor, Kjk the kinetic and Wjk the gravitational po- 
tential tensor. The dots in (|4"B")) indicate other possible terms due to magnetic 
field, pressure etc. If an external force is exerted on the object, we have to add 
to the right hand side of equation (jig)) the term [52, p. 280] 

Vjk = -~ / p\ xk^r— + Xj^r— d*x , (49) 



2 J \ dxj dxk 

where <£> e xt is the external potential and the case of a cosmological constant 
corresponds to 

$cxt = -^Ar 2 . (50) 

Therefore the new virial theorem which accounts for the cosmological con- 
stant [53] takes the form 

^=4K jk +2W jk + h jk A+.... (51) 

This is, in principle, a differential equation for Ij k provided Kj k and Wjk are 
given or corresponding differential equations in these variables are derived. It is 
very often more convenient to consider a less demanding task by simply noting 
that the trace W of Wjk is negative whereas the trace K of Kjk is positive defi- 
nite. Then the gravitational equilibrium i.e. d Ij k /dt = yields the inequality 

~AI+\W\>0, (52) 
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where I denotes the trace of the inertial tensor Ijk- To appreciate the meaning 
of this inequality we specialize to the case of constant density p = p. It is then 
easy to show that ([52"!) takes the form 

8717? > AA , (53) 

where the quantity A depends only on the geometry of the object under con- 
sideration. More specifically it reads 

167r f r 2 d 3 x . . 

J p 

where $n is the Newtonian part of the non-relativistic gravitational potential. 
For spherically symmetric objects one easily calculates A = 2 and therefore the 
virial inequality is simply 

47rp>A. (55) 

We already mentioned the similarity between Buchdahl's inequalities (which are 
a consequence of hydrostatic equilibrium) and the inequalities describing the 
validity of the Newtonian approximation. A second kind of such a coincidence 
(if one can call it coincidence) occurs also here with regard to equation (|55|) . 

and (fT4")l . To add to this coincidence we mention that the authors of [32] by 
demanding stability of circular orbits derived further the inequality 

471-p > 4tt/9 c > A . (56) 

No doubt, all three inequalities originating from different premises have some- 
thing in common. The similarity between the three results is not trivial. Cer- 
tainly, the hydrostatic and gravitational equilibrium are intuitively related. 
However, the hydrostatic equilibrium necessarily involves pressure whereas (|55|) 
is derived without using this concept. Indeed, the virial equation for an object 
in gravitational equilibrium with pressure reads 

\kl + W + 2K + ZPV + ... = 0, (57) 

o 

where the pressure P is constant over the volume V (note that P can have 
different signs depending whether the pressure is internal or external) and the 
dots indicate the presence of terms related to magnetic fields, rotational effects, 
etc. We think that this issue deserves a closer inspection. 

With the help of the virial theorem we can show a drastic astrophysical 
effect of the cosmological constant. Inequalities like (1521 and (|53"|) are useful to 
check if a given astrophysical system is in gravitational equilibrium without the 
knowledge of the kinetic tensor or, which is the same, the average velocity (v 2 ) 
of the components of the system. It is often of utmost interest to know this 
velocity. The knowledge of the shape of the object, its density profile and A 
allows one to calculate 



\W\_ _ 87T Pvac 

M 3 M 



I. (58) 
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To appreciate the effect of A let us assume a constant density and the shape of 
the astrophysical object to be an ellipsoid. Then the W and I can be calculated 
analytically [52] . For the inertial integral / we get 

I = P^a^a-z (aj + a\ + af) , (59) 
where a, are the axes of the ellipsoid. The mean velocity can be now written as 

(v 2 )ciii P soid = ^-pPvacaia 2 a 3 (a\ + a\ + a§) ^^r e iii pso id - 1^ . (60) 
For the oblate case (ai = a 2 > a 3 , e = y/l — a\/a\ ) one calculates 

Toblato = 7^ , 2 • (61) 

The prolate case (oi = a 2 < a 3 , e = \/l — a\/a%) gives 

1 prolate — 2 ' V Uz V 

! + 2 fe) 6 



The triaxial case (a\ > a 2 > »3, = cos 1 (a 3 /ai), k — y/ (a\ — a|)/(af — a|)) 
is simply 

_ a 1 a 2 F{9,k) 

triaxial — o , _2 , _2 ' ^1 

where fc) is the elliptic integral of the first kind. Note now that for a 
flattened prolate ellipsoid we can approximate 



prolate 



\ 3 / 11 

ai \ ( I + e 



Since the nowadays preferred value of p V ac is (0.7 — 0.8)p C rit we can say that 
if the constant p/p C rit is, say, 10 3 , it suffices for the ellipsoid to have the ratio 
a\/az ~ 10 in order that the mean velocity of its components approaches zero. 
This is valid always under the assumption that the object is in gravitational 
equilibrium. This effect is due to the relatively large cosmological constant. 
In general we can say that in flattened astrophysical systems in gravitational 
equilibrium, the mean velocity gets affected by the cosmological constant. The 
denser the system, the bigger should be the deviation from spherical symmetry 
to have a sizable effect. It is interesting if such an effect can be observed in 
reality which would confirm the existence of A. One can paraphrase this also by 
looking at the results from a different perspective. If we are certain that a given 
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astrophysical object is in gravitational equilibrium, then equation (52) and (53) 
would put a stringent bound on A in the case of strong deviation from spherical 
symmetry. 

From the last equations we can infer that the power of equilibria concepts 
with A lies in analyzing relatively low density astrophysical objects whose shapes 
deviate from the spherical symmetry. It is evident that this can easily be done 
using the virial theorem. We have shown that for spherically symmetric objects 
the virial theorem and the hydrostatic equilibrium yield surprisingly similar 
results. It would be therefore desirable to have a tool to probe the hydrostatic 
equilibrium with A for arbitrarily shaped bodies as well. The full examination 
of this aspect of hydrostatic equilibrium goes beyond the scope of the present 
paper. We can, however, outline a direction in which such an investigation 
might proceed following [54] . The basic equations for the Newtonian concept of 
hydrostatic equilibrium are 

VP = -pV* , V 2 $ = 4np - A , (65) 

together with (|26[) as the equation of state. Defining u = p/po, £ = sx with 
s 2 = 47T/9 2_7 /k7 we obtain for u = u(£) 

(!-7)(-^) +VI\nu = a A u 1 -^ -u 2 -> , (66) 

which is a generalized Lane-Emden equation. For A = we have to set oa = 
whereas the case A ^ requires cla = 2. Furthermore for A = 0, po is arbitrary, 
but is fixed to be po — p vajc for non-zero cosmological constant. We will examine 
the consequences of equation (|66p elsewhere. However, here we can already 
mention a difference between the cases A = and A ^ due to the fact that 
the cosmological constant sets a scale for the density. Consider the case 7 = 1. 
Then with 7 = 1 and for A — we have a scaling property in the following 
sense: if u(£) is a solution of the equation, so is A 2 m(A£) [55]. This similarity 
property is lost if a a = 2. 

Coming back to the virial theorem we note that interestingly one can elimi- 
nate the cosmological constant if we address the question whether an astrophys- 
ical object is in gravitational equilibrium. Let the system which is in gravita- 
tional equilibrium be described by I\, Wi and K±, etc. This system serves us as 
a reference system. Defining 621 = I2/I1 = IA12 and provided A ^ 0, a second 
system (denoted by the subscript 2) is also in gravitational equilibrium if the 
following relation is fulfilled 

( ' 21 - W 1 +2K 1 +E 1 > (6?) 

where 3j contains terms related to pressure, magnetic field, etc. On the other 
hand, if we know two systems bounded by an external medium at a constant 
pressure P cx t (this means that 3j = — 3P eK tVi assuming that there is no magnetic 



15 



field) are in equilibrium, we can solve for the pressure and obtain 



P. 



V 2 



W2 
Wi 



£21 



2K-i 



K2 



£21 



(68) 



Note that the salient and necessary assumption for equations (IBTl) and is a 
non-zero cosmological constant A. With these equations we can put in relation 
two systems in equilibrium or solve for the variable common in both systems. 
Solving for the kinetic energy of system 2 we have 



1 



K 2 = - (2e 2 iK 1 



£21 Wl -W 2 + £21^1 - S 2 ) • 

Since this is a positive definite quantity, we find an inequality for W 2 

1 



\W 2 \ > e 2 i [\Wi\- 2K t - Si t- — =2 

£21 



(69) 



(70) 



which gives an estimate of the left-hand side. 

As an example, we assume again that the system 1 is in steady state and 
that both systems are in equilibrium with an external pressure P e xti &s before. 
Under these conditions we can solve for the pressure and find 



< 



\W- 



w 2 



3Fi \\Wt\ 



£21 



V 2 

£21 

Vx 



(71) 



which can serve as a lower bound for the external pressure. The question that 
rises here is, what can we identify as the reference system? On the one hand, one 
could argue that at every different astrophysical scale there must be a standard 
system, so that one must find, for instance, a standard galaxy in order to obtain 
information about a certain other galaxy. The exploration of these interesting 
aspects will be the subject of further research. We conclude this section by 
the remark that the virial theorem with non-zero A has also been successfully 
applied in cosmology [56]. 



8 Scales of the cosmological constant 

To appreciate the orders of magnitude estimate we reinstall here the Newtonian 
coupling constant G which we put to 1 in the preceding sections. The cos- 
mological constant sets scales for density, length, mass and time. The density 
scale in equation (1) indicates the so-called coincidence problem. The latter 
can be formulated as the question as to why we are living in an epoch in which 
Pvac ~ Pcrit- More precisely, since p cr it is epoch dependent, the coincidence 
Pvac ~ Pcrit is a time coincidence of a sharp transition from small ratio Q\ to 
fl\ ~ 1 [57]. We can reformulate this coincidence in terms of length, noting 
that 

r A = —= = 0.476 H^ 1 = 1.43 x 10 3 h^ ^pc , (72) 
v A 
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which is of cosmological order of magnitude and close to the radius of the visible 
universe. The combination 



1 



= 3.61 x 10 22 /in 1 M, 



Pcrit 



(73) 



M A 



Pvac 



is close to the mass of the universe. As noted in [31] this is not trivial since A 
as a constant could be either smaller or bigger than what is assumed today. We 
could equally well live in an epoch where the universe is smaller or bigger or 
the mass of the universe is not close to M\. Interestingly, both such scales play 
an important role in the considerations of the validity of the weak field approx- 
imation and in the hydrostatic equilibrium in the form of Buchdahl inequalities 
(see the discussion in section [4|. Hereby the differences of the numerical factors 
are of order one. Had we lived in a different universe or a different epoch of 
our universe, the above scales would appear only in the Newtonian limit and 
the Buchdahl inequalities, which by itself is a remarkable technical coincidence. 
However, in our universe these scales appear in all considerations: cosmology, 
hydrostatic equilibrium and Newtonian limit. Notice also that one of the equa- 
tions of motion of a test body in the Schwarzschild metric with A is the relation 
between the affine parameter A and the time t which reads 



with rs — GM . The case E = sets the limit of validity of the coordinate 
system. This means that the expression on the left hand side of (|74|) is zero 
which in turn leads to the solutions r m ; n = 2r$ and r max = a/3Ya- Albeit the 
cases r — r m i n , r max are singular points due to the choice of coordinates, the 
condition above tells us also that we should not go beyond the horizon of a 
black hole and also not beyond the horizon of the universe. The latter emerges 
only because of the above coincidence of length mentioned above. In the same 
context, i.e., considering the motion in a Schwarzschild spacetime the scale R 
defined by 



not only yields astrophysically relevant length scales (this is due to the combi- 
nation of a small and large quantity), but also has the meaning of the largest 
possible radius within which bound orbits in the Schwarzschild metric are pos- 
sible [58]. 

Curiously, not all scales are of cosmological or astrophysical relevance even 
though they have the right order of magnitude. The time scale connected to A 
is 



which is close to the age of our solar system, but should be considered as an 
accidental coincidence. 




(74) 



R 3 = 3r s r 2 A 



(75) 



T A = 4.65V 1 Gyr , 



(76) 
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The scale = y/X = 3.1 x 10 42 GeV has no connection to particle physics, 
but interestingly establishes one more of Dirac's large numbers in the form 
m p /(GVA) - 3.2 x 10 41 where m p is the nuclcon's mass. 

9 Summary 

In this work we investigated the astrophysical limits on the cosmological con- 
stant for spherically symmetric bodies. These limits emerged from examinations 
of hydrostatic and gravitational equilibrium as well as from constraints on the 
validity of the Newtonian limit. Although the resulting inequalities are based 
on quite different premises, they do resemble each other to an extent which 
let us conjecture a deeper underlying principle or connection. As a side result 
we conjecture the existence of a new class of static and spherically symmet- 
ric perfect fluid solutions with the cosmological constant. By coincidence the 
length and mass scales set by A appear as the horizon and mass of the universe 
and simultaneously as the limits of validity of the weak field expansion and the 
hydrostatic equilibrium. This adds to the puzzle of the coincidence problem. 

In this context we note that the inequalities based on the hydrostatic equi- 
librium were derived in the general relativistic framework whereas the virial 
inequalities required only a Newtonian approximation. We suspect that gen- 
eral relativistic virial equations might be a more suitable tool to compare the 
hydrostatic equilibrium with the virial one. 

Already the Newtonian virial theorem is a powerful equation in many sit- 
uations where the astrophysical object is not spherically symmetric. Indeed, 
flattened, large and diluted objects cannot be in gravitational equilibrium if the 
cosmological constant is as large as claimed nowadays. It is therefore of some 
interest to extend the consideration of the hydrostatic equilibrium allowing for 
an arbitrary shape of the objects. Unfortunately the general relativistic treat- 
ment of arbitrarily shaped objects is quite out of reach at the moment, but we 
have indicated how to do that in the Newtonian approach. This will be covered 
elsewhere. Similarly some of our results, like the inequality (43) still await its 
practical exploitation. 

The inequality 2p > p vac which we found from various considerations is valid 
provided that the object is spherically shaped and in equilibrium. Assuming this 
and taking into account that p can be at the most of the order of the background 
density i.e. p cr it, this is not a weak inequality. Indeed, it comes very close to 
the preferred p vac ~ 0.7p crit . 

In this work we have laid down the basics of equilibria concepts including a 
cosmological constant. As indicated in the paper one should continue this study 
by examining hydrostatic equilibrium for non-spherical objects and generalize 
the equilibria concepts also to other theories of the Dark Energy [59] . 

We can look upon many inequalities derived in the main text by imposing 
the condition of equilibrium from two different angles. Assuming a known value 
of A the inequalities limit the possible values of the density of the astrophysical 
object if the latter is assumed in equilibrium. If we know that an object is 
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in equilibrium and assuming we can infer the mean density of such an object 
from some other considerations, then the same inequalities put a limit on the 
cosmological constant. 
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